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Dependance of weak localization cone with the photons brownian Doppler shift 
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We report the first observation of the dependence of the coherent backscattering (CBS) enhanced 
cone with the frequency of the backscattered photon. The experiment is performed on a diffusing 
liquid suspension and the Doppler broadening of light is induced by the brownian motion of the 
scatterers. Heterodyne detection on a CCD camera is used to measure the complex field (i.e. the 
hologram) of the light that is backscattered at a given frequency. The analysis of the holograms 
yield the frequency and the propagation direction of the backscattered photons. We observe that 
the angular CBS cone becomes more narrow in the tail of the brownian spectrum. The experimental 
results are in good agreement with a simple theoretical model. 
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Coherent backscattering (CBS) of light is a photon self 
interference effect which leads to an enhanced intensity 
cone in the backscattering direction. This effect, which 
is related to the solid state physics phenomenon of weak 
Anderson localization, has been extensively studied both 
experimentally and theoretically. Since the first obser- 
vation of CBS from colloidal suspension [l], Q , the phe- 
nomenon has been observed on many experimental sys- 
tems such as powders [3!,|3,[5[, liquid crysta ls M , photonic 
crystals [7,], and cold atoms gases [1 i, [MJllThe CBS 
effect has also been studied on acoustic [l^, [iS 13 and 
seismic waves |l5l]. 

In most of the CBS experiments done in the optical do- 
main, the angular distribution of the backscattered light 
is measured by using an incoherent detection method. 
The signal that is detected by a photomultiplier, a pho- 
todiode or a CCD camera is proportional to the optical 
field intensity. In this letter we propose to study the 
CBS effect more precisely by using a coherent detection 
method [16|. The backscattered field is summed with a 
coherent reference field in order to detect the interference 
pattern of the two fields. The detection is thus sensitive 
to the phase of the backscattered field. 

In this letter, we propose a detection scheme for the 
CBS effect, which involves a heterodyne detection of the 
light on a 2D detector (CCD camera). Since the pixel to 
pixel relative phase of the detected backscattered field is 
accuratel y in easured, we record the digital hologram of 
the field [l3l- By adjusting the frequency of the hetero- 
dyne detection local oscillator, we are able to precisely 
select the frequency of the backscattered field that is de- 
tected. We are thus able to measure the backscattered 
light average intensity as a function of both the propaga- 
tion direction (as done commonly) and the frequency i.e. 
as a function of both the wave vector direction k/|k| and 
the modulus Ikl. The CBS effect can then be observed 
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FIG. 1: Experimental Setup. L: main laser, MM: moving 
mirror, M: mirror, LO and SIG: local oscillator and signal 
arms, BE: beam expander, AT: attenuator, AOMl, AOM2: 
acousto optic modulators, , BS: beam splitter, S: cell filled 
with the diffusing suspension. A: rectangular aperture, CCD: 
CCD camera. 



with a backscattered field whose frequency is shifted with 
respect to the incident incoming field. Since this shift is 
a Doppler shift which is related to the motion of the scat- 
terers, we can study the influence of the motion of the 
scatterers on the CBS signal. We will see in particular 
that the angular width (and more generally the angular 
shape) of the enhanced backscattered cone depends on 
the frequency offset. 

Fig [T] shows the experimental setup. The 780 nm, 50 
mW beam of the laser L (Sanyo DL-7140-201, current: 
95 mA) is split in two beams LO and SIG, which are 
both vertically linearly polarized. The main beam SIG 
(> 90% of the power, frequency z^l) illuminates the dif- 
fusing sample S (7 cm high x8 cm width xl cm depth 
cell which is filled with a pure or a diluted intralipid 
suspension: Endolipide 20%; B. Braun Medical SA, dif- 
fusion anisotropy factor g = 0.65 at A = 780nTO 18]). 
To detect the S backscattered light signal, the LO beam 
(< 10% of the power) is mixed to the signal beam with 
the beam splitter BS. The LO + SIG interference pattern 
(or heterodyne beat) is recorded by a CCD camera (PCO 
Pixelfly: 12 bits, vccd = 12.5 Hz exposure time 20 ms 
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FIG. 2: A) Linear scale 1024 x 1024 calculated image of the 
average k-space field intensity \E\^ for ulo = i^l- One pixel 
is 1.13 X 10~* radian. B) Cut of the image along the dashed 
line (black curve). Angular gain of the experiment (solid grey 
line). 



with 1280 X 1024 pixels of 6.7 x 6.7 fim) and transferred 
to a PC computer. 

To avoid saturation of the CCD camera the LO beam 
is attenuated by the attenuator AT. The power is finely 
adjusted by moving the mirror MM. By using the AOMl 
and A0M2 acousto optic modulators (Crystal Technol- 
ogy: i'AOMi,2 — 80MHz), the LO beam frequency vlo 
can be freely adjusted: vlo = vl + (vaomi - vaom2)- 
Two numerical synthesizers with the same 50 MHz quartz 
reference clock provide the i'aomi,2 radiofrequency (RF) 
signals. Two beam expanders BE (x25: Spectra Physics 
model 334) enlarge the LO and SIG beams. The BE fo- 
cus are adjusted to get a plane wave LO beam on the 
CCD, and a plane SIG illumination beam on the sample. 

To get the complex field E at frequency u, we use the 
4-phase heterodyne variant of the phase-shifting method 
[ITI Il9j . We record sequences of 4 CCD images li at 
ij = to-'ri/vccD with i = 0..3. E interferes with the local 
oscillator field Elo, and /, = l-Bioe^^'"'^"*' -K£'e-''2^'^*' P 
where p — —1. We can choose v^o ^ v — vccd/^, and 
\Elo I ^ 1^1 • E \s, then given by the simple 4-phase equa- 
tion E — a[(/o— /2)+j(/i— /a)] where a is a constant. Due 
to the motion of the scatterers, the backscattered light 
frequency spectrum is broad, and the 4-phase equation 
performs the heterodyne detection on the two sidebands 
frequencies u^ = vlo ± '^cc£>/4, yielding to the field 
components E{v^). The 4-phase equation should be thus 
rewritten as: E{v+) + E*{v-) = (/q - h) + j-ih ~ h) 
[20J . In the holographic reconstruction, the 2 compo- 
nents E{i'^) and E*{v~') yield to 2 different images cor- 
responding to the holographic grating orders ±1, which 
may spatially overlap. To spatially separate the images, 
we have (i) reduced the pupil size by using a black rect- 
angular aperture A (1 cm width x5 cm high), (ii) moved 
the sample at a respective distance Dsampie — 25 cm to 
the CCD, and (iii) shifted the LO beam off axis {9 on 

FigO m ■ 

We have recorded N = 256 series of 4 images (4 x A^ 
images) with the 1 cm thick cell filled with Intralipid 



10% (optical depth ~ 20 Is* where Is* is the transport 
mean free path) . The LO beam is unshifted in frequency 
(i.e. i^Lo = vl), and the CCD acquisition time is 20 
ms per image. For each series, we have calculated the 
real space field E{x,y) by the 4-phases equation, and 
the k-space field by 2D Fast Fourier Transform (FFT): 
E{kx,ky) = FFT E{x,y). The angular distribution of 
the field intensity is then obtained by summing \E\'^ over 
the N series of images. Figl^K) shows the k-space field 
calculated average intensity \E\'^. Calculations are per- 
formed on a 1024 x 1024 matrices obtained from the 
1280 X 1024 CCD raw data after truncation. From now 
on, we will consider the angular size ao of the k-space 
pixels as the units for the angle. Due to the FFT we 
have: ao = A/(1024 x dptxei) = 1-13 x lO^'' radian. 

The two vertical oblong spots (1,2) seen on Fig[21^ are 
the two grating orders (±1) images of the sample. Due 
to the off-axis configuration, images (1) and (2) do not 
overlap. The images are blurred (the images are oblong 
while the aperture A is rectangular) since the holographic 
reconstruction (FFT) is done at infinite distance while 
the object is located at a finite distance Dsampie- The 
(1,2) spots also represent the angular distribution of the 
backscattered light, and we can see narrow brighter spots 
(3,4) in the middle of the spots (1,2) that are the CBS 
enhanced cones. Due to the aperture A and moreover 
the fact that the sample is quite far away from the CCD, 
the light, which is backscattered far from the illumina- 
tion axis, does not reach the CCD. This means that the 
angular detection efficiency of our setup is not flat (the 
wide spots 1 and 2 do not fill the k-space). To analyze 
the shape of the CBS cone, we have made a cut along the 
vertical line that passes in the center (xg, yo) of the CBS 
white spot (dashed line of Figl2] A) and we have plotted 
the CBS signal (black curve of Figl2]B). To improve the 
SNR (signal to noise ratio), the plotted curve is obtained 
by summing over the lines a; = xq — 5 to xo + 5. The 
CBS peak (3) is clearly seen, but, as mentioned above, 
the background is not flat. 

We have corrected the signal from the background dis- 
tortion by using the measured holographic data to calcu- 
late the complex field image of the sanrple in the aper- 
ture A plane. From that calculation we got the position 
of the aperture A i.e. Dsampie^ the shape of the illunn- 
nation zone, and the local variations of the illumination 
intensity. To get holographic signal \E\'^ that we should 
obtain if there is no CBS effect (i.e. if the photons are 
scattered by the sample in all directions with equal prob- 
ability), we have calculated the blurred intensity image 
of the objet that is expected with a focusing at infinite 
distance. We have by this way the expected background 
signal i.e. the efficiency of our setup for detecting the 
backscattered light in a given angular direction. This 
quantity is plotted as a grey solid line on Figl2j3, which 
agrees with the observed shape of the CBS background. 

The heterodyne detection technique allows the study of 
new physical effects since the frequency of the backscat- 
tered light that is detected can be finely adjusted within 
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FIG. 3: Angular shapes of the CBS peak for \vlo ~ vl\ = 
0(a), 100(b), 200(c), 300(d), 400(e), 500(f) and 600 Hz(g). 
The curves are obtained by summing over the x = a;o — 5 to 
a;o + 5 lines. To make all curves visible on the same diagram, 
logarithmic scale is used on the y vertical axis {y = 0.05 to 
1.0). Horizontal axis is angle (1.13 x 10" radian per point). 
A) Experimental results. B) Theoretical predictions. 
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FIG. 4: Frequency spectrum of the backscattered light (log- 
arithmic scale vertical axis). Points are the CBS background 
measured on curves a to g of Fig|3j\. Curve 1 (solid line) is the 
brownian theoretical shape calculated by Monte Carlo, curve 
2 (dashed line) results from [22] ), curve 3 (solid grey line) is 
the single scatterer Lorentzian. Horizontal axis is vlo — vl- 



the brownian spectrum. We have repeated the previous 
experiment for different LO frequencies: u^o ^vl = +0, 
+100 ...+600 Hz (curve a, b, ...g of FiglSlA.). To correct 
the angular detection efficiency of our setup, the mea- 
sured signal (solid line of Figl2j3) is divided by the ex- 
pected background signal (heavy grey line) , which is cal- 
culated from the vlo — vl experimental data, in order 
to gain the highest SNR (signal to noise ratio) measure- 
ment. Fig. [31A shows the CBS cone angular shape ob- 
tained by this way. The detection efficiency correction 
is seen to be efficient and the background is flat. As ex- 
pected, the CBS cone signal and the background both 
decrease when when |i^LO ~ vi\ increases (i.e. when the 
detection is done in the tail of the brownian spectrum). 
Moreover, one observes a new physical effect, the angu- 
lar width of the CBS cone also decreases with \vho — Vh\ 
(curve a to curve g). 

This effect can be understood quite simply. By chang- 
ing the frequency of the detection (i.e. Vho\ one changes 
the distribution of travel paths that contribute to the sig- 
nal. Since each scattering event broadens the photon fre- 
quency spectrum, one selects, by pushing vj^o in the tail 
of the brownian spectrum, travel paths with many scat- 
tering events (to get photons far from the center of the 
brownian spectrum, one needs more scattering events), 
i.e. longer travel paths. This corresponds to a narrowing 
of the CBS peak. 

To make a quantitative analysis of this effect, we have 
measured the brownian diffusion constant Db of the in- 
tralipid particules by considering that the CBS back- 
ground measured on FiglSj^ yields the shape of the brow- 
nian frequency spectrum. In the case of the backscatter- 
ing by single diffusor, this spectrum is a Lorentzian of 
width Ak^Ds [2311 ■ In the multi diffusion case, each travel 



path yields a Lorentzian of width DB^i\f\i^ (where 
Si is the sommation over the scattering events i, and 
q.j = ki — ki_i is, for each event, the incident versus exit- 
ing difference of wavevectors [22]). The spectrum is then 
a sum of Lorentzians. By considering a set of 10000 travel 
paths, we have calculated the shape of the multi diffusion 
spectrum (FigUcurve 1). The paths are obtained by sim- 
ple (scalar approximation, Henyey-Greenstein scattering 
function) Monte Carlo simulation [2J| with g = 0.65. The 
multi diffusion spectrum (curve 1) has been compared 
with the single scatterer spectrum (curve 3), and with the 
diffusion wave spectroscopy predictions (curve 2). This 
one is the Fourier transform oig\(t) = eT'^'^^^ -^^'^ [22| 
with 7 = 1.9 (i.e. for parallel polarisation and g — 0.65). 
We have adjusted the frequency scale to make curves 1 
to fit the measured point (a to g). The width of curve 3 
(^k'^DB = 74 Hz) yields Db- 

We have used our Monte Carlo set of paths to calcu- 
late the angular shape of the CBS cone. For each path, 
the photons may travel forward (from the first to the last 
scatterer) or in reverse. The interference of the two paths 
yield a cosine angular intensity distribution analog to a 
Young fringe pattern. By summing these fringes over 
the paths, one gets the CBS enhanced cone. This intu- 
itive idea dates to the CBS original papers [l|, 0, and 
is illustrated by a web Java applet [2a|. We have gen- 
eralized this calculation by considering the angular and 
frequency distribution of light for each path, and then 
summing over the paths. The CBS curves of FigO cor- 
respond to the 2D cuts (angle, frequency) distribution 
along the vj^o — i^l — 0, 100, ...600 Hz lines. We have 
plotted on Figl3j3 the curves that result from the ran- 
dom walk calculation. To make a good theory versus 
experiment comparison, the Fig|3}3 curves are obtained 
by summing over the 11 cuts made on the experimental 
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FIG. 5: Same as Fig|3] but with arbitrary units linear scale 
vertical axis. The points are experimental results, the solid 
grey curves the theoretical predictions. 



data (x = lo — 5 to xo + 5). The angular position of 
the CBS cone, the width of the CBS peak, and the CBS 
enhanced factor rj have been determined from the curve 
a of Fi^3J^ i'^LO — vl). The other curves (b to g) are 
then obtained without any free parameters. 

As seen, the experimental results (Figl3J\_) are in good 
agreement with the theoretical predictions (FigjSjB). The 
CBS enhancement factor [rj = 1.48) is lower than 2 be- 
cause the linear polarizations are parallel, and because 
we have summed over 11 parallel lines. With good opti- 
cal component and careful matching of the BE focus, our 
system could be diffraction limited by the CCD area. In 
that case, the angular resolution should exactly be one 
pixel of the k-space. We have plotted the CBS curves 
with a linear scale for the vertical axis. For each curve 
(a to g) the scale is arbitrarily adjusted. 

In this letter, we report the first observation of a de- 
pendence of the coherent backscattering (CBS) enhanced 
cone with the frequency of the backscattered photons. 
The narrowing of the CBS cone that is observed in the 
tail of the brownian spectrum is in good agreement with 
a simple model which describes the effect of the travel 
path distribution. We do not observe any break of reci- 
procity because the scatterers do not move during the 
travel of the photon and the scattering cross section does 
not depend on the motion of the scatterers (as it could do 
in cold atom experiments [li|). These results are made 
possible by the development of a new detection method 
of the backscattered photons, which is based on hetero- 
dyne detection of the photons with a CCD camera. This 
method, which allows the study of both the angular and 
spectral properties of the backscattered light, opens the 
way to the observation of new physical effects in disor- 
dered systems. The authors thank D.Delande, A. Tourin, 
J. de Rosny, G. Montambaux and D. Bonn for fruitful 
discussions. Authors acknowledge also french ANR for 
its support. 
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